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I. INTRODUCTION 



The interpretation of neutrino-nucleus scattering data relies on accurate knowledge of 
the vA cross section. Above the quasielastic peak higher resonances increasingly contribute 
to the cross section. Neutrino experiments generally measure integrated yields which in- 
clude contributions from quasielastic nucleon knock-out as well as from higher resonance 
production. We discuss, in the following, the effect of A production in neutrino scattering. 
It has already been shown in the case of electron scattering that a qualitative description 
of the data can only be achieved when A-h excitations are included [P],EJ. Charged-current 
neutrino-nucleus scattering has been used as tool to investigate axial-vector form factor of 
the nucleon ||. For low momentum transfers the Q-dependence of the form factor can be 
parameterized by a dipole mass Ma, Eq. ( |A19|) . We will consider the effect of A excita- 
tions on the experimental extraction of Ma- Kim et al. || have examined nuclear structure 
corrections to the extraction of Ma but do not consider the A. Singh and Oset have 
included A-h but calculated in nonrelativistic formalism. 

The article is organized as follows. First, in Section |I| we introduce the relativistic 
mean-field formalism of the nucleus including the A resonance, deriving the appropriate set 
of nuclear response functions. In Section [TT1] we discuss the results for cross sections and 
yields as modified by the inclusion of the A channel. 

II. FORMALISM 

In this section, we derive the inclusive cross section for quasielastic charged-current 
neutrino scattering including A-h excitations in the nucleus. We consider a neutrino with 
four- momentum k=(E u , k) which scatters from a nucleus via W ± boson exchange producing 
a charged lepton with four- momentum k' '=(-Ek' M') ■ Using an impulse approximation and a 
Fermi gas description of the nucleus the formula for the double differential scattering cross 
section for mass number A is given by (we assume a symmetric N = Z nucleus): 



2 



d 3 a AG 2 F cos 2 # c |k'| 

rf^ k ,^ k , " 32 **pE v \ n u A ) , (1) 

where p = 2k F /3ir 2 is the baryon density with Fermi momentum kp. Gp denotes the Fermi 
constant and 9 C is the Cabibbo mixing angle. The leptonic tensor is defined as 

Lfj,u — 8 (jt^k u + k v k^ k ■ k g^ u =p itaPfivk k ^ , (2) 

with the minus (plus) sign denoting neutrino (anti-neutrino) scattering. U A U is the polariza- 
tion tensor of the target nucleus for the charged weak current. Here we consider p-h, rig, 
and A-h, n A ^, contributions to the polarization: 

nr = ng + nft • (3) 

The expressions for the p-h polarizations have been derived in previous publications [§]||. 
The weak interaction contains vector current (v) and axial-vector current (a) contributions. 
Therefore we split n A ^ into: 

nft = (W + W + W + (W ■ (4) 

(n™J'" / and (n^) M ^ are interference terms of the vector and axial-vector currents. In a 
Hartree approximation the polarization tensor can be written in the form 

(n A ,V = -*E/ 4^nn,(-^-p) s " a (p) r UviP) G(p- q ) ] 



n,p 



+ (q l M^-q t t) (i,j) = (a,v). (5) 

We rewrite the interference term as 

(nftr = (n^r = «r a0 q a n% h . (6) 



S^ip) is the Rarita-Schwinger form of the free spin 3/2 propagator with momentum p 



p 1 + M A 



- M A + ie 



1 „ „ p*y -p v ^ 



ST - - n~T~FT + 



(7) 



3 ' ' 3 Mi 3M A 
Note that this expression is not unique and other forms of the A propagator have been 
considered (see Ref. 0). The differences enter the off-shell behavior of the propagator 



which do not affect the following calculations. The vector part of the nucleon-delta vertex 
has been studied in the case of the 7iVA transition ffl, 



r(q,p)=V2F A T ± 



(8) 



The isospin raising (lowering) operator originating from W + (W ) exchange is defined 
through 

T ± = ^=(T 1 ± l T 2 ) (9) 

where 2x4 isospin matrices T l satisfy || 

T\T ] ) j = 5 ij - ~rV . (10) 
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The form factor Fa is defined in the appendix [Eq. ([A 1 8| ) ] . 



The vertex for the axial iVA transition is given by |||9| 

T ± 

T % = , (11) 

with the axial form factor Ga, Eq. ( |A19| ). The parameter r^A indicates the strength of 
axial iVA transition and will be discussed later. In the noninteracting limit the nucleon 
propagator G(p) reduces to the free fermion propagator G°(p) for a relativistic Fermi gas 
with Fermi momentum kp. We consider only the density-dependent part G° F as vacuum 
contributions do not enter at the Hartree level. In the rest frame of the nucleus one obtains 

G° F (p) = (jf+ M)—5(po - E p )9(k F - |p|) . (12) 
F P 

Using an impulse approximation, the imaginary parts of (II^/J enter the cross section 
Eq. ([J). As long as the A is assumed to be stable, the imaginary parts can be calculated 
analytically. The resulting expressions are given in the appendix. 

In a relativistic mean field description of the nucleus, nucleons and As interact with the 
background of scalar (a) and vector (u) meson mean fields. The interactions are assumed 
to be analogous to those of nucleons with possible new couplings, g s A and g A . The values of 



the couplings can be constrained somewhat by fitting the quasielastic peak of A production 
in electron scattering (1].|2|] . 

In a relativistic mean- field approximation (MFA), the noninter acting nucleon propagator 
G° F {p) is replaced by 

Grip) = (f + M*)—6(p - E p )9(k F - |p|) (13) 
E p 

where 



M* = M - S N , E* = Jp 2 + M* 2 , p*" = j/* - V N g^° . (14) 



The scalar (Sn) and the vector (Vjv) self-energies can be obtained for a given fcp JTO 
Analogously, the mass and the momentum in the A propagator S^ v are replaced by [|[] 

.\/ A - .\/ A • .S A . f -> i*" =e- . (15) 

The calculation proceeds in the same way as in the case of the free A. The expressions for 
the imaginary parts are given in the appendix. For simplicity we assume the A self-energies 
5a, Va to be the same as nucleon self-energies SV, Vat- 

In free space, a A decays into 7rN with a width V = 115 MeV. In the medium the 
situation is more complicated. The 7rN decay channel is partially suppressed because of 
Pauli blocking, i.e., the phase space available to the nucleon produced in the A decay is 
reduced by the Fermi sea. However, the A in the medium has additional channels of decay 
and obtains a "spreading" width with the main decay mechanism from A + iV — > N + N. 
These two competing effects cancel each other partially. In our calculation, lacking better 
theoretical and experimental knowledge of the A width in nuclear matter, we assume a value 
identical to the free width T = 115 MeV 



As we are mainly interested in integrated cross sections, the results do not depend 
strongly on the way the width of the A is treated. We adopt a simple method to include 
the decay width by averaging the nuclear response over the A mass with a Breit-Wigner 
distribution The averaged cross section follows as 



< 



d 3 a 




roc 

(/,)/(m 2 )/ / >W) 

JM 2 



(16) 



k' 



m 2 d 2 Q ]i >dEi L f 

M A r 



(Mi - /x 2 ) 2 + m 2 r 2 



(17) 



integrating from threshold to infinity. 



III. RESULTS 



In this section we present result for the A-h calculations of charged- current neutrino 
interactions. We discuss the case of muon neutrinos but the general features of the results 
hold for electron neutrinos as well. In addition to a relativistic Fermi gas calculation we 
consider the effects of the mean field and include the decay width of the A. The target 
nucleus is assumed to be 16 with a Fermi momentum kp = 225 MeV. For the strength of the 
axial NA transition a simple argument using the A decay width suggests a value r Nl \ ~ 2.2 
|| whereas constituent quark models give a somewhat smaller value of r^A = ~ 1-7. 
We choose an intermediate value r^A = 2 for our numerical calculations. 

Figure [I] shows the double differential cross section for measuring an outgoing muon 
produced by an incoming neutrino with energy E u = 1 GeV and three-momentum transfer 
|q| = 0.5 GeV. First, note that the curve neglecting the decay width of the A has a peak 
around go = 0.37 GeV which agrees with the expected elastic A peak at 



assuming the initial nucleon at rest. The A cross section is similar to the p-h cross section. 
Therefore, measurements of integrated quantities cannot neglect A production. Including a 
finite delta width reduces the peak height by about 30 to 35 percent but does not significantly 
reduce the total integrated strength. 

Nuclear matter effects are included using a mean- field approximation (MFA). Here we use 
the same scalar and vector couplings for nucleon and A (known as "universal couplings"). 
For k F = 225 MeV, a self-consistent nuclear- matter calculation yields the effective masses 




(18) 
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Ml = 931 MeV, M* = 638 MeV, and the vector self-energy V N = V A = 239 MeV. Mean- 
field results are also shown in Fig. [1|. p-h and A-h contributions are reduced by about 30 
percent. Both peaks are shifted to higher energies due to the smaller effective masses. 

As we have seen that A-h excitations can contribute substantially to the charged-current 
cross sections, it is interesting to study the effect of A production in neutrino-nucleus scat- 
tering experiments. Experiments using muon neutrino beams have measured 

f(E u ) denotes the spectrum of the neutrino beam. da/dQ 2 is given by 

do [ Q < 7T d?a 

dQ 2 Jo E u \k'\ dEvd 2 Q k , % 1 ' 

where Q 2 = —q 2 = q 2 — g 2 and the cut off for the energy transfer Q c reads 

q 2 -m 2 E v m 2 , 
AE V q 2 — 

Using the neutrino spectrum from the charged-current experiment at BNL the resulting 
cross sections are shown in Fig. |2| where we used an axial mass M4 = 1.09 GeV. At larger 
momentum transfers the contribution from A-h excitations is as large as the nucleon knock- 
out. 

We now discuss the influence of deltas on the extraction of the axial mass from quasielas- 
tic data. Different kinds of experiments are possible. If only a charged lepton is detected, 
all A events will be included on an equal footing with p-h excitations. Alternatively, an 
experiment could detect pions and thereby separate A events producing real pions from p-h 
excitations. However, a significant fraction of A excitations lead to two-particle two-hole 
excitations without a real pion. A A in a nucleus can decay via A + iV — > N + N. This is 
related to either pion absorption or weak meson exchange currents (involving an intermedi- 
ate A). It may be difficult to separate two-particle two-hole from one-particle and one- hole 
final states. The axial mass is often fit to reproduce the Q 2 dependence of observed events. 
This cancels some errors from unknown flux normalizations. Therefore it is interesting to 



consider the Q 2 dependence of (a) p-h excitations only (b) p-h plus that fraction (see below) 
of A-h excitations leads to 2p-2h (c) p-h plus all A-h excitations. An incorrect value of Ma 
could be extracted if one assumes only p-h excitations while the data is "contaminated" by 
significant A-h excitations. 

We leave it to the experimental groups to analyze their data in detail. For example, 
Ref. |§ analyzed their data assuming only p-h excitations and extracted a value of Ma = 1-09 
GeV with a very small statistical error of ±30 MeV. To estimate the uncertainty in this 
extracted value of Ma from A-h excitations we try and fit the Q 2 dependence from 0.3 - 1 
GeV 2 of our full calculation (p-h plus some fraction of A-h) with a p-h only model. Thus 
the full calculation assumes some value of Ma (which is essentially arbitrary) and the p-h 
calculation attempts to reproduce this result by using a possibly different value of Ma- The 
important quantity is the difference between the assumed and extracted Ma- This may 
represent some of the systematic error (from A-h excitations) in a p-h only analysis of data. 

Figure |3| (a) shows a likelihood function for reproducing our theoretical results assuming a 
p-h only free Fermi gas with different values of Ma- (Note, all of the theoretical calculations 
used Ma = 1-09 GeV.) For theoretical calculations assuming only a p-h response (solid 
line), the input Ma = 1-09 GeV is of course reproduced in the fit. However, for theoretical 
calculations including either all of the A-h (dashed line) or half of the A-h events (dots), 
Ma is underestimated by 70 to 90 MeV. 

This factor of half represents a very crude estimate of the A + iV— > N + N to A — > Nir 
and A + iV — > N + N branching ratio. Theoretical results [|T2| , p~3[1 are consistent with this 
factor. However, there could be both important Q 2 and model dependence in this branching 
ratio. Further theoretical work on the branching ratio would be very useful. Alternatively 
one could try and measure it in coincidence electro-excitation experiments. 

Finally in figure |3| (b), we fit theoretical calculations including scalar and vector mean 



fields (assumed independent of momentum) as described in Ref. [14J] with a free p-h cal- 
culation without mean fields. Again large shifts in the extracted Ma are found. A events 
tend to increase the effective cross section at high Q 2 which might be fit with a smaller Ma- 
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Alternatively, mean field effects tend to reduce the cross section at high Q 2 which can be 
fit using a larger Ma- Thus there is some cancellation between the two effects. However 
this cancellation is unlikely to be perfect and the theoretical uncertainties are large. We 
conclude that the theoretical uncertainty on an extracted Ma could be of order 0.1 GeV and 
thus large compared to the claimed experimental error of ±0.03(stat)±0.02(syst) GeV 0. 

IV. SUMMARY AND OUTLOOK 

We have calculated charged- current neutrino cross sections including A-h excitations of 
the target nucleus. The calculation was done for free deltas as well as including the effects 
of relativistic scalar and vector mean fields in the nucleus. A-h excitations are found to give 
significant corrections to quasielastic nucleon knock-out processes in experiments measuring 
neutrinos in the GeV range. In extracting the axial form factor of the nucleon from neutrino 
scattering data the A-h channel enters with similar strength as p-h contributions. This may 
introduce significant error in the extracted nucleon axial form factor. 
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APPENDIX: 



The imaginary parts of the A-h polarizations are evaluated analytically and listed below. 
Note that the second term of the polarization, Eq. (^), with (gy, — ► — q a ) vanishes for on-shell 
A: 

Im(n^)^ = ^^IIV + V(M 2 + Ml) - (M 2 - M\ ) 2 ]£ A , (Al) 



Im(n- ) 00 = -j^rjFl 



187r|q 



4g 2 i?3 A + 4go/3^ 



3 ^ 

4 o„2/j\ / f2 i j/2 \ i / 7\ /f 2 7i/f2\2 , j j\ #2 J2l jjA 



+ [q 4 - 2q\M 2 + M£) + (M - M^) + 4ilTg 2 ]£ 1 z 
37r|q 

2a n2 .E* + 2q Ef + q*E* 



(A2) 



Im(n^)^ = ^G^f , (A3) 



a 



97r|q| MJ 



A 

Im(n^) 01 = -_^_G^[2gb^ + (3gg + q 2 + M 2 - M 2 )£ 2 A + g (/? + 2q 2 )£ A ] , (A5) 

Im(n^) 11 = --^—(^Et + Qo((3 + 2q 2 )£ A + + g 2 + M 2 )£ A ] , (A6) 
Im(nX h ) = --^G A F A [2q 2 E^ + q (3E*] , (AT) 

where a — [(M + M A ) 2 - q 2 ] and (3 = q 2 + M 2 - M\. Also 

= ~ — (n = 1,2,3), (A8) 

n 

£ A = Min(£ F , E Amax ) , (A9) 

/3g + |q|V/3_ 2 -4MV 
2g 2 



^Amm — ^ • l AiU J 



Im(n^) 22 is obtained from the relation, 

99 n 00 -n u -rF u 
n 22 = & . (ah) 

In the mean field approximation, the A and nucleon masses in the propagators are shifted 
by strong scalar fields. Since we take the same interaction as in free space, the polarizations 
involve complicated traces. After a little algebra, the polarizations are written as 

8 r E F r 1 In I T^ v 

(nSZr = -^/ dE P d x f\ ( - V2 vv m2m . , (A12) 

3 Jm* J-i 8it z [p + q) 1 — M A + le 
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8 r E F r 1 Inl T^ v 

8 r E F r 1 Inl if 111 "* a T 

(DSr - / jf , ^ p * (p + ; )2 ^£ +i6 , (AM) 

where T^, and T m result from evaluating the traces. These can be determined straight- 
forwardly as 

?T = — (p*-t*-MM*) 



3M A 



- t ■ qM A M* A { q n* u + q v t*") 
+ (2M* A t ■ q - M A t* ■ q){ty v + pV)M* 
+ + t u t*^)(M* A M A q 2 - 2t ■ qt* ■ q) 

-2{t ■ qfig^Mf - t*H* v ) + 2t»f[(t* ■ qf - q 2 M* A 2 } 
+ {qH™ + gT") M A [M A (M* A 2 p* ■ q + 2t* ■ p*t* -q)+t- qM* A 2 M] 
+ M A M A [{t»q v +t V) (t* -qM-p* ■ qM* A ) + (q^p* u 
+q u p*^)(t* ■ qM A — t ■ qM* A )\ 
+ M* A M A [{m* v + t u t*^)(p* -qt* -q- MM* A q 2 ) 

+{p*H* v + p* u t*^)(t -qt* -q- q 2 M A M* A )] 
+ 2(g^q 2 - q^q u )M 2 A M* A 2 (t* • p* + M* A M) 
+ 2t-qp- qM A M* A {g» v M A 2 - p*»p* v ) 
- 2t* ■ p*Ml(g^( t* ■ q) 2 + q 2 t*H* v ) 

+ 2M A M* A t* ■ qg^ v (t* ■ p*t ■ q - M A M* A p* ■ q - 2t ■ qMM* A ) 
+ M A M* A (q 2 M A 2 - (t* ■ q) 2 Wp* u + fp*")) , (A15) 



TZ = [3(t* 2 - M* A 2 ){p*H* u + p* v t^) - 16t*H* u (p* ■ t* + M*M* A ) 

6M A 

+g ^t* 2 (p* ■ t* - 2M*M* A ) + 15# M V • t*M A + 18g flu M*M A 3 } , (A16) 



T va = -77^(E* - ^X)[2p • qM* A - M* 2 M A - 4M*M* A M A + q 2 M A 
6M A |q| 

r*2 



3M A M A ] . (A17) 
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The vector form factor of the NA vertex is given by 

-< M * + M > 9 'l ^VYl ^V" 2 (A18) 



M((M A + M) 2 - g 2 ) 2 V 0.71GeVV V 3.5GeV 



which has dipole form including some phenomenological corrections |15|]. For the axial vertex 
we use the nucleon axial form factor 

1 26 
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FIGURES 

FIG. 1. Double differential cross section for (v, /x) scattering from 16 showing p-h and A-h 
contributions separately. Results are shown neglecting (solid line) and including the width of the 
delta (dashed line). The thin lines denote cross sections in the mean-field approximation. 

FIG. 2. da/dQ 2 averaged over BNL antineutrino spectrum are shown for Fermi gas (solid line) 
and mean field approximation (dashed line), p-h and A-h results are shown separately. 

FIG. 3. Likelihood function (normalized) of fitting experimental antineutrino-scattering results 
varying the extracted axial mass, assuming a free particle-hole (p-h) only Fermi gas response. Parts 
(a) and (b) show results for Fermi-gas and mean-field approximations, respectively. Curves are 
shown for p-h events (solid line), the sum of p-h and A-h events (dashed line), and p-h plus half 
of the A-h events (dots). Note, all theoretical calculations assumed Ma = 1-09 GeV. 
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